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^ ■ Abstract : In this paper we study several means of compensating for thermal lensing 

(N : 

5_i . which, otherwise, should be a source of concern for future upgrades of interferometric detec- 

■ tors of gravitational waves. The methods we develop are based on the principle of heating 
00 ' 

^ ■ the cold parts of the mirrors. We find that thermal compensation can help a lot but can not 

I . do miracles. It seems finally that the best strategy for future upgrades ( "advanced configu- 
> ; 

^ . rations") is maybe to use thermal compensation together with another substrate materials 

^ . than Silica, for example Sapphire. 

o : 

^ , 

o ■ 

1j ; I. INTRODUCTION 

U '. 

Large scale interferometric detectors of gravitational waves such as LIGO [011 or VIRGO 
'. I@l currently under achievement. Their first data taking should begin in a few years 

I from now. But from their planned sensitivity and from what is currently known from 

astrophysical sources, direct detections of gravitational waves is far from sure to happen 
with this first generation of interferometers. Plans for improving the sensitivity are already 
being developped, and the design of second generation interferometers is well advanced 
0. Among the noises that have to be lowered, the shot noise is one of the fundamental 
ones. It is indeed the dominant one in the high frequency region (typically above 1 kHz) 
of the sensitivty of first generation interferometers. The shot noise improvement does not 
amount to naively increase the power of the laser source for example, since the mirrors of the 
interferometer have low but non zero absorptions. This results in thermal gradients in the 
mirrors and a possible degradation of the detector performance. Thermal effects in mirrors 
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have been shown to be under control for the present design of interferometric detectors 
but they are hkely to be one source of concern for future upgrades and even to 
introduce a new source of noise via the couphng of laser power fluctuations to absorption 
asymmetry in the interferometer arms . Indeed, lowering the shot noise limited sensitivity 
by one order of magnitude means increasing the circulating powers by two. The planned 
circulating powers in VIRGO for example are of the order of 15 kW stored in the kilometric 
cavities and about 1 kW crossing the (Silica) substrates of the cavity input mirrors. The 
measured absorptions in the Silica susbtrates are about 1 ppm/cm for the Heraeus 311SV 
P], and the best absorptions in the coatings have been measured to be about 1 ppm pHlT 



With these figures, there will be about 25 mW dissipated in the input mirrors (10 mW due 
to absorption in the susbtrate and 15 mW due to absorption in the coating). This gives 
reasonable losses (mainly due to the thermal lensing effect) which does not affect too much 
the sensitivity. But, it is an other story if the circulating powers increase by as low as a 
single factor 10, for example if we keep the interferometer as it is, and simply upgrade the 
laser to a 10 times more powerful one (upgrade from 10 W to 100 W in the case of VIRGO). 
In the current design of "Advanced LIGO" |Q, a laser source of about 100 W is planned 
and hundreds of kW are planned to be stored on the high finesse arm cavities (about 500 
kW in case of Silica mirrors up to about 800 kW in case of Sapphire mirrors). Clearly, at 
a time when R&D is going to be funded for advanced configurations of gravitational-wave 
detectors, means have to be investigated in order to control the induced thermal lensing 
which otherwise could be desastrous. 

The thermal lensing takes its source in the apparation of temperature gradients in the 
mirrors, these gradients being due to the non uniform beam intensity crossing the mirrors 
or hiting the coatings. In this article, we develop some means to decrease these gradients. 
The basic idea is very simple: we have to heat the "cold" parts of the mirrors in order 
to homogenize the temperature field. This can be done by an external light source with 
a proper profile, at a proper wavelength, and being absorbed by the parts of the mirror 



2 



that are not heated by the main (about Gaussian) beam itself, or by a thermostat system 
that fixes the external temperature of the mirrors to some value (principle of the 'electric 
blanket'). Note that the implemantation of such a concept is planned for "Advanced LIGO" 



r^. The goal of the article is of course not to provide a technical solution to this problem 



but rather to indicate the prospects of the idea of compensating for temperature gradients 
by heating the mirrors. First we will recall how we estimate the temperature field and the 
induced thermal lens for a mirror heated by a Gaussian beam either by absorption in its 
substrate or by absorption in its coating and suspended in vacuo. We will then give some 
orders of magnitude of the lensing effect and related losses for first generation interferometer 
and then for an upgrade. We will take the VIRGO numbers for the numerical evaluations, 
but without really a loss of generality. For the upgrade we will assume simply a 100 W laser 
instead of a 10 W laser so all the circulating powers have to be multiplied by 10, so 10 kW 
crossing the input mirrors and 150 kW stored inside the long cavities. In a second part, 
we will consider the case of a mirror heated by an extra beam with a ring-shaped profile 
as the one roughly produced by a battery of diodes located all around the mirror (radial 
heating ogf the mirror). We will then consider the case of a beam coaxial with the main 
beam and again with a ring-shaped profile (in this case we have to take into account the 
damping of this beam due a priori to a strong absorption along the mirror axis). Finally 
we will consider the case of a mirror with its lateral side "in contact" with a thermostat 
of fixed temperature. In each case we will compute the losses of the mirror, for a one way 
propagation of a perfect TEMqo beam through the mirror (and so experiencing the thermal 
lensing), for various absorbed powers from the external light source or various thermostat 
temperatures (the main parameters of the studies), in order to look for an optimum (minimal 
losses). We will also consider briefiy the case of Sapphire substrates, instead of Silica ones, 
but the substrate materials choice will not change our conclusion. We will finally discuss 
the infiuence of the (main) beam size. 
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II. TEMPERATURE AND THERMAL LENS IN A MIRROR 



A. The mirror, the Heat equation and the boundary conditions 

We consider an long cavity input mirror: a Silica substrate of width /i ~ 10 cm and 
radius a ~ 20 cm with a reflective coating located in the intra-cavity side. The mirror is 
heated by a 'laser' beam, assumed in this part to be a perfect TEMqo beam with a waist 
u^o — 2 cm either in the substrate, where the circulating power is Psub, or in the coating, 
where the power is Pcoat- The coating, modelled as a very thin layer, is located at z = —h/2 
(see fig.l). 

We will note a the lineic absorption coefficient in the substrate and e the absorption 
coefficient in the coating. If alone in vacuo, the mirror can only loose heat by radiation. 
Fortunately the temperature rise in the mirror is supposed to be moderate so the radiative 
heat flux can be linearized (this will be assumed in all the paper): 

F = (7(r^ - T^) ~ AaT^(T - T^) (2.1) 

where Tg is the temperature of the vacuum tank (supposed uniform) and a is the Stefan 
constant (corrected for Silica emissivity). In this part, starting from now, we will note T 
not for the absolute temperature, but its deviation from Te (T — Te — ^ T); it is valid since 
everything is now linear. The steady-state Heat equation to solve for the mirror is then: 

- KAT = al{r) (2.2) 

together with the boundary equations: 

dT 

-K—(r = a,z)=AeaTo^T(r = a,z) (2.3) 
or 

dT 

-fC— (r, z = h/2) = AeaTg^T{r, z = h/2) (2.4) 
dT 

-K—{r, z = -h/2) = el{r) - AeaTj^T{r, z = -h/2) (2.5) 
oz 

The linearisation of thermal equations (boundary conditions) implies that we can study 
separatly the different sources of heating: absorption by the coating or the substrate of the 



Gaussian beam or of any other beam. Note that, in Eq. (|2.2|) , we have neglected the possible 
attenuation of the beam through the mirror. This is of course valid since we consider here 
only the main laser beam for which the substrate is a low loss one. We will have to take 
into account this effect for auxiliary beams if their wavelength is within absorption bands of 
the Silica. 



B. Absorption in the coating 

If we consider only absorption in the coating, then the Heat equation (|2.2|) becomes the 
homogeneous one AT = 0, a general solution of which being expressed as a Dini series: 

T(r, z) = Yl (^me'^'"' + Bme-^-') Ukmr). (2.6) 

m 

which respects well the cylindrical symmetry of the problem. The first boundary condition 
(2.3) shows then that the coefficients C,m = ci ^ are the zeros of the function 

xJi(a;) - rJo(a;) = 0, (2.7) 

where r is the reduced radiation constant: r = AeaT^a/K. The functions Jo{kmr) are then 
a complete set of orthogonal functions, for functions defined in [0,a], with scalar products 
given by |15| 



h{kmr)3o{knr)rdr = ^1 + Jo(Cm)^'^nm- (2.8) 

The coefficients Am and B^, are then given by the two last boundary conditions (2.4) and 
(2.5), and we find finally the temperature: 

r(r .1 - V ^P-^p-'^^VSa (Cm-r)e-^"^'^-^^ + (Cm + r)e-^"'^ 

where the pm are the coefficients of the Dini expansion of the intensity /(r) 

I{r) = T.PmMkmr) (2.10) 
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The Prn are then given by, taking account the normahsation condition of Eg. ( |2.8|) : 

2 1 1 

Pm = — rFT^TTTT^ / Hr)h{kmr)rdr. (2.11) 



In this section we consider only the heating by the main beam (assumed to be a perfect 
TEMoo of power P) 

I^r) = ^e-^'"/< (2.12) 

TTW^ 

SO that Pm can be computed as 

r) = — I I Q-klw^o/^ (2 13) 



at the condition a >> wq 16 



The optical path distortion due to thermal lensing is then given by |T4 



f+h/2 

^{r) = — Tir,z)dz (2.14) 

ai J-h/2 

where dn/dT is the thermal index coefficient of the substrate materials. After straightfor- 
ward integration, we find 



= §^Ef^ ,, ^ I , -U H ^oikmT). (2.15) 

dT K ^ km (Cm + ^) - (Cm " ""^^ 



Some quantitative results and figures are given below. 



C. Absorption in the substrate 



In this case we have to solve the full Eq.( p.2|) . For this purpose, we add to the general 
solution Eq. ( |2.6| ) of the homogeneous equation a particular solution. A particular solution 
of Eq. (pl^) is easily found to be J2m "^^Jo(^m^) so that the complete solution is 



T{r,z) 



^ (Ame'-^ + Bme-'-^ + ^) Jo{kmr) (2.16) 



The boundary condition (2.3) gives again that x o, is the m-th zero of xJi(x) — tJo(x) 
and the two others (in fact the same) give the coefficient Am = Bm- The temperature is 
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finally 



nrt ( (p-k,n{h/2-z) , -km{h/2+z) \ 



By using Eg .( p. 141) , we find also the thermal lens 



dn ^ apm / , 2r 1 — e \ / , n / x 



D. Numerical examples with VIRGO and VIRGO 'upgrade' 

We give in this part some examples of typical temperature fields and thermal lensing in 
the VIRGO context. As the thermal theory is linear, we first give results separately for the 
cases of absorption in the coating and in the substrate for unity absorbed power. In order 
to quantify the impact of thermal lensing on the VIRGO performance we compute, as a 
figure of merit, the losses £ induced by the thermal lens. These losses are computed as the 
deviation from unity of the coupling coefficient of a perfect normalized (power unity) TEMqo 
beam propagating through the mirror, so experiencing the thermal lens, to the same TEMqo 
beam. The coupling coefficient is nothing but the squared modulus of the inner product 
between the two beam complex amphtudes. We have then 

r°° 2i7r'P(r) 2 

/ ^00(^^)6 ^ iJoo{r)27irdr 



(2.19) 



where '?/'oo(^) = \/2/7rwQexp(— r^/ti^g) is the amplitude of the normalized TEMqo beam. 

On figure 2 we can see the temperature fields and the corresponding thermal lens profiles 
for 1 W absorbed respectively in the coating and in the substrate. We recall that, although 
the temperature fields are different for the two cases, with a maximal temperature rise 
around 10 K per absorbed Watt in the coating and less than 3 K per absorbed Watt in 
the substrate, the thermal lenses are very similar, with, in both case a maximal optical 
path difference (OPD) about 3/im per absorbed Watt. In the planned case for VIRGO, we 
will have ePcav — 10^^ x 15 x 10^ ~ 15 mW absorbed in the cavity mirror coatings and 
ahPsub — 10^^ X 10 X 10'^ ^ 10 mW absorbed in their substrates, so the corresponding 



maximal OPD is about 7.5 x 10"^ /im. If, in case of upgrade, for example the laser power 
is increased by one order of magnitude, all the absorbed powers increase accordingly. It will 
be then rather 100 mW that will be absorbed in the coating and the same amount in the 
substrate, so the corresponding maximal OPD is about 0.75 /im. If now we compute the 
induced losses in the case of VIRGO, we find £ ~ 1.7 x 10~^, which is reasonable, and in the 
case of the upgrade C ~ 0.16 which is disastrous. The situation seems a little better if we 
discard the parabolic part (obtained by a fit of the thermal lens in the mirror central part) 
from the thermal lens. Indeed, a parabolic lens effect can be compensated for by a proper 
matching lens at the input of the interferometer. Of course, this works only if the thermal 
lensing is well symmetric between the input mirrors of the kilometric cavities, otherwise it 
will be impossible to compensate for parabolic parts of both thermal lenses by some input 
matching optics. In this case, the losses become C ~ 2.5 x lO^'' in the case of VIRGO and, 
in the case of VIRGO upgrade, £ ~ 2.2 x 10~^ which is still too high. These results clearly 
call for decreasing the thermal gradients in case of (even minimal) upgrade. 

NB: in the following, all the numerical results involves the "VIRGO upgrade" case {Psub — 
10 kW and Pcoat — 150 kW) and the up-to-date values for the absorption coefficients (e ~ 1 
ppm and a ~ 1 ppm/cm), except the final ones given for the "advanced LIGO" configuration. 

III. COMPENSATION OF THERMAL EFFECTS BY ABSORPTION IN THE 



In this section, we will try to solve the thermal gradient problem by heating the external 
part of the mirror by a ring-shaped beam, in addition of the main Gaussian beam, which is 
always supposed to be absorbed in the coating and in the substrate. We model the intensity 
of the extra beam as 



SUBSTRATE OF A RING-SHAPED BEAM 




if r < Oe 




(3.1) 
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where /q is constant. The total power in the beam is then 

Pe = /ovr(a2 - al) (3.2) 

If the ring-shaped heating zone is obtained with the help of a a number of diodes located all 
around the mirror, then the quantity a — represents nothing but the penetration length 
(related to the hneic absorption of the diode wavelength) of these beams in the substrate. 

There is thus an extra source of heating in the substrate with respect to the previous case. 
The temperature field due to absorption of light in the substrate is then the superposition 
of the one computed in section 2 and the one due to absorption of the ring-shaped beam. 
For the later, the solution is essentially the same, except that a has to be replaced by the 
lineic absorption for the wavelength of the extra beam, and the coefficients of the Dini 
expansion of the intensity have to be changed. Let's note p^-* the coefficients related to 



/e(r). According to Eq.( |2.11| ), we have then 



By direct integration of x^q{x) in x^i{x) and by substitution of /q by the corresponding 
power Pg; we have finally 

(e) -^e '^Cm 1 f i f \ 



= .ia--al)r^ + ClHCmY V'^^^^ ~ -^^Ma) ] (3.4) 



From Eq. (|2.17|) , we can then write the total temperature field due to absorption in the 
substrate: 

run -I- A- 'n(^) / (p-k„^{h/2-z) , -k„^{h/2+z) \ 

Tir, .) ^ E (l - rL-^-^^±l-^) MK,r> (3.5) 

In the same way, we obtain the thermal lensing profile by substituting apm by apm + ctePm 
in Eq.( p.l8D . The expressions for the temperature field and thermal lens due to absorption 



in the coating are of course unchanged. 

The impact of the supplementary heating by the extra beam depends on two parameters, 
the radius which indicates which size of the mirror is heated and the power of the beam. 



or, better, the absorbed power in the mirror aehPg. In order to optimize the situation, 
figure 3 shows the mirror losses as a function of the radius a^, the absorbed power being 
fixed for each value of to the value that minimizes thes losses. We see clearly that there 
is an optimum at — 2 cm, that is the numerical value of the waist wq of the main beam. 
Henceforth, we set then Og = 2 cm. The figure 4 shows the mirror losses as a function of 
the absorbed power. We see then that the heating by the ring-shaped beam dramatically 
improves the situation: the losses can be decreased from the (inacceptable) value £ ~ 0.16 
down to £ ~ 2.7 x 10~^, so a 'gain' of about a factor 60 for the losses. An extra factor about 
2 (corresponding to losses as low as 1.2 x 10~^) can be gained if we substract the parabolic 
part of the thermal lens. Heating the cold parts of the mirror seems then of great benefit, but 
the obvious drawback is that a large amount of light power has to be dissipated, about 38 W 
for the optimal case (that gives a dissipated intensity /q — 3 W/cm^), that means a source 
of yet higher power. If we imagine a battery of laser diodes around the mirror, delivering say 
1 W each, at least 38 of them will be needed. In order to understand why so large powers 
are needed, let's take a look at the shapes at the profiles of the thermal lenses in the optimal 
case and for too low or too high absorbed powers (see fig. 5) We note that the total thermal 
lens has a quite perfect fiat profile at the center of the mirror. This compensation for the 
thermal lensing due to absorption of the main beam can be only achieved for a particular 
power absorbed from the extra beam. If the absorbed power is too low or too high, then 
there remains a concave or convex lens. 



IV. COMPENSATION BY ABSORPTION IN THE SUBSTRATE OF AN 
ATTENUATED RING-SHAPED BEAM 

We now consider the physically different case of an external ring-shaped beam propagat- 
ing parallely to the main beam and being (strongly) absorbed by the mirror substrate. We 
just saw in the previous section that, in order to be efficient, the external beam must deposit 
a large amount of power in the substrate. That indicates that in the case now considered 
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the extra beam will be strongly attenuated along its propagation through the mirror. We 
have then to solve the heat equation with an extra source of heat being the absorbed light 
intensity. As in the previous section we adopt a ring-shaped profile for the intensity but 
with now an attenuation along z. We suppose that the new light source is in the intra-cavity 
side, but it can be in the other side, it will give exactly the same new thermal lens. The 
intensity can be then written as 

{0 if r < ai 

(4.1) 
I^Q-»2(z+h/2) \{ ai<r <a2 

where a2 is the lineic absorption of Silica at the wavelength of the new light beam. Note 
that by chosing wavelengths near 0.5, 1.4 or 3 /im (OH absorption bands), we can have 
values of a2 in the range 0.1 to 0.95 /cm [[l^,!!!] We note also that we can have a priori an 



external radius 02 of the beam lesser than the mirror radius. 
We have now to solve the Heat equation: 



K^T = a2le{r,z) (4.2) 



The general solution of the homogeneous equation has always the form of Eq. ( |2.6[ ). We have 
to add a particular solution of ( |4.2| ) to the later. Such as a particular solution can be looked 
for under the form I] tme~'^^^Jo(^m'") and we find: 

= ?T/^e--V^ (4.3) 

The solution is finally: 

T(r, z)=J2 (Ame'-' + B^e-'-^ + ^-/^^^e-'^^^^^'+A Jo(A;„r). (4.4) 
The boundary condition Eq.(2.3) g fcm X 0, is the m-th zero of xJi(x) — 
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rJo(x), and the two others Eq.(2.4) (with e = 0) and Eq.(2.5) give: 



^ «2 _,,V2 («2Q + r)(r + Cn^) + (^^2^ - r)(r - C^)e-("^+^"')" 

Kkl-al"" (^^ + ^)2_(^^_^)2e-2/c™/^ 



The thermal lens due this temperature field is derived from Eg. ( |2.14]) : 



(4.5) 



(4.6) 



dT ^ K{ki- ai) (Cm + ^) - (Cm - r)e ^ ^ ^ ^/ 

(4 

We obtain then the total temperature field by superposition of the field due to absorption 
of the main beam (in the coating and in the substrate) and of the field due to absorption 
of the auxiliary ring-shaped beam. The same for the thermal lens. Let's notice that the 
absorbed power in the substrate from the auxiliary beam is 

Pabs = J «2e-"^("/'+^)Pot^^ = Po{l- e-"^") = a2hPo i \ (4.8) 

where Pq = 2n J I^rdr = /ovr(a2 — af). 

There are now four free parameters in this part : the extra beam inner radius ai and 
outer radius 02, the power of the source Pq and the absorption 02 at the wavelength of the 
source. The two last ones can be merged into a single one, that is the absorbed power. We 
will see indeed that the results actually do not depend on Pq and 02 but simply on Pats- 

As a first example, figure 6 shows the losses induced by the total thermal lens as a 
function of Pabs for cti = 2 cm, 02 = a and 02 = 0.2 /cm, that is a penetration factor 
1 — exp(— ^2^) — 0.865. We note that the minimal losses are (again) of the order of 2.7 x 10^^ 
and are obtained with Pats — 38 W. If we allow now a2 to vary, we find in fact the same 
kind of results : minimal losses about 2.7 x 10"^ for the same absorbed power about 38 W. 
If we optimize now with respect to ai and 02, we find, just as in the previous section, an 
optimum at ai ~ 2 cm (corresponding to the size Wo of the main beam) and at 02 — 19 cm, 
just below the mirror radius a ~ 20 cm. 

12 



The minimal losses achievable by this method are L ^ 2.6 x 10~^, and are given by 
ai ~ 2 cm, ai ~ 19 cm and whatever the value of ai- Of course the higher ^2, the smaller 
the needed source power Pq. It has to be noticed that the parabolic part removal from the 
thermal lensing doesn't improve the minimal losses here. There is actually a gain when the 
absorbed power is well below or above, but not around the optimal value (about 38 W), for 
which the lens is almost perfectly fiat at the center of the mirror and parabolic removal can 
not improve anything anymore. The minimal losses quoted above are then an absolute limit 
for this method. 

The conclusions of this section are then similar to the previous one, except that there can 
not be any extra gain due to parabolic lens removal. In both case, in the optimal absorption 
case, the mean temperature of the mirror is about 20 K. The approximation of low heating 
(T ^ Te) is still reasonable (error on the radiative flux calculation less than 10%). 



V. COMPENSATION BY FIXING THE MIRROR CIRCUMFERENCE 

TEMPERATURE 

A. Heat Equation general solution and boundary conditions 

We suppose in this section that the temperature of the mirror circumference is fixed to 
a temperature To by some thermostat (for instance an electric blanket or some equivalent 
device). We always suppose that the vacuum tank is at the temperature Tg, that is the 
temperature 'seen' by the mirror sides at 2; = ±/i/2. From now on, we will take Tq as the 
reference temperature, and T(r, z) will then denote the deviation from this reference. The 



heat fiux,Eq. (^TT| ) is in this case: 

F = a ((To + Tf - Tf) ~ + ^aTl(T^ - T,) (5.1) 

We will then have the same expression for the radiative flux on the mirror sides, except 
that appears a supplementary (constant) term that we will note C = 4aT^{TQ — Te). The 
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boundary conditionds can now be expressed as: 
dT 

-K—(r = a,z) = (5.2) 
or 

dT 

-K—(r, z = h/2) = AeoTj^Tir, z = h/2) + C (5.3) 
oz 

dT 

-K— (r, z = -h/2) = d{r) - AeaTj^T{r, z = -h/2) - C (5.4) 
We solve the Heat Equation as in section 2. A general solution of the homogeneous Heat 



Equation is always given by Eq.(2^). The first boundary condition Eq.(5.2) implies that 
the km are now defined as (m = km x a are the zeros of Jo(a;)- The functions Jo{kmr) form 
again a complete orthogonal set with the normalization condition given by [|r^ 

/ Jo{kmr)Jo(knr)rdr = -a^Ji(Cm)^^nm- (5.5) 
JO 2 

The Gaussian beam intensity can always be expanded on this set, and the expression of the 
Pm has to be modified in accord the new normalization constants : 

7ra^Ji(Cm)^ \ 8 J 
with the same notations as in section 2. 

In the two other boundary conditions, the constant C can be expanded on the set of the 

Jo{kmr). If we note C = J^Cm^oikmr), we have: 

2 



a^Ji(Cm) 



^2 







CJo{kmr)rdr. (5.7) 



We note that this amounts to compute the Hankel transform of a radial window function, 
and after a simple calculation, using that the derivative of xJi{x) is xJo(x), we find 



Cm Jl (Cm) 

In practice, the number of terms to consider in the Dini Series will be fixed by the 
condition 

l^j2^7^jf\^re[0,a[. (5.9) 

From figure 7, we see there is an optimum at ~ 700, and, consequently we will chose 
A^ = 700 in the following. Note that for all the previous numerical results, we have chosen 
A^ = 50, which was actually comfortable. 
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B. Absorption in the coating 



In the case of absorption in the coating, the boundary condtions Eq.(5.3) and Eq.(5.4) 
become: 

-K (k^A^e^-"-'^ - k^B^e-''-^/^) = ^ (A^c'^-'^/^ ^ B^e-'^-'^l^) + (5.10) 

— K(h A p-'^^VS _ L R pkmh/2\ _ _ ll^ ( A , pk^h/2\ _ 



(5.11) 



The solution of the system is then, after a straightforward calculation. 



n _ -k,„h/2 i^Pm ~ Cm) (Cm + ''")— CmjCm ~ ^""^ /p. , o\ 

" (Cm + ry - (Cm - r)^e-2^™'^ ^''''^ 

So we have the solution of the Heat Equation for heating in the coating by absorption of 
the gaussian beam by the Dini series Eq. (|2.6|) . 

The thermal lens is again calculated from Eq. ( p.l4D and from the temperature Dini 
expansion: 

Of course, if C = 0, that is all the vanish, we find again the expression of Eq.( |2.15| ), the 
only difference being then the definition of the km- 



C. Absorption in the substrate 

The solution of the Heat Equation Eq. (|2.2|) is again the superposition of the general 
solution of the homogeneous equation AT = and a particular solution. The solution 
has again the form of Eq.( p.l6D , and, this time, the two boundary conditions Eq.(5.3) and 



Eq.(5.4) give: 



rapm _|_ acm 

(Cm + r) - (Cm - r)e-^™'^ ^ ' 
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so that the temperature is: 

' apm 

m \Kkl (Cm + r) - (Cm -r)e 



rapm _|_ ac„ 

(5.16) 



After integration of Eg. ( |2.14| ), we find the thermal lens profile: 



(5.17) 



\ Kkl^ km (Cm + t) - (Cm " r)e 
If all the Cm = 0, we find again the same expression as Eg. (p.l8| 



D. Numerical results 

The only free parameter in this compensation method is the temperature Tq of the 
thermostat (if the temperature Tg of the surrounding vacuum tank - the tower - is fixed). 
The figure 8 shows the total losses as a function of the thermostat temperature. We note 
there is a minimum occuring for Tq ~ 311 K, so a temperature rise about 15 K with respect 
to the tower temperature. But the minimum is rather high, about £ ~ 1.2 x 10^^, so one 
order of magnitude larger than the minimal losses obtained with the other compensation 
methods described in previous sections. Clearly this method is not competitive. Removing 
the parabolic part does not help much here and changing the temperature Tg does not change 
roughly the results. For sake of completeness, figure 9 shows the temperature and thermal 
lens profile for the optimal thermostat temperature (Tq ~ 311 K). We note an average 
temperature mirror less than 10 K and a roughly (but not sufficiently) flat profile of the 
thermal lens at the center of the mirror 



VI. SAPPHIRE SUBSTRATES 

It is already well known that replacing Silica substrates by Sapphire ones can already 
improve the situation. For example, for the VIRGO 'upgrade' situation, the losses due to the 
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thermal lensing can be decreased from £ ~ 0.16 to £ ~ 1.2 x 10~^ if we simply replace Silica 
by Sapphire. This is simply due to the much larger conductivity of Sapphire {K ~ 33.0 
W.m^^.K^^ vs K ^ 1.38 W.m"^.K^^ for Silica). If we use in addition a compensation 
method, we can have even better results. In the case of compensation by absorption of a 
ring-shaped beam (sections 3 and 4), the losses can be decreased down to £ ~ 1.7 x 10^^, 
that is 17 ppm (1.1 x 10^^ - 11 ppm - with parabolic part removal). The drawback is that 
this optimum is reached with a larger absorbed power Pabs — 46 W than in the Silica case 
{Pabs — 38), but it is not a dramatical increase of the needed extra-source light power. 
Note that in the case of compensation with the help of a thermostat (section 5), we can 
achieve losses less than 10"'^, but this requires a fine tuning of the thermostat temperature 
by a fraction of degree in excess from the tower temperature, something non realistic (for 
example if compared to the tower temperature inhomogeneities) . 

VII. INFLUENCE OF THE BEAM SIZE 

Increasing the beam size is also a way to decrease the intensity absorbed in the mirrors, 
and so a way to decrease the thermal gradients. In this section, we are going to study the 
influence of the (main) beam size on the thermal losses and, in particular, the impact of a 
compensation method. Note that in the "advanced LIGO" concept, large beam sizes are 
planned to be used, typically wq ~ 5.5 cm in case of Silica substrates |Jl2|]. We show only 
results for the compensation method by absorption of a ring-shaped beam (see sections 3 
and 4), since the last method (section 5) is not competitive. For each beam size, the minimal 
losses are found (1) by optimizing the inner radius of the ring-shaped beam (2) by optimizing 
the absorbed power. We find that the minimal losses are roughly independant of the main 
beam size, but, what is changing is the optimized absorbed power from the extra beam. The 
minimal losses are £ ~ 2 x 10"^ (if wq > 3 cm) and the needed absorbed power decreases 
down to about 4.5 W for a beam size wq = 6 cm, instead of 38 W for a beam size wo = 2 cm, 
as shown by figure 10. Increasing the beam size then does not change the minimal possible 
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losses, but it changes the needed power to reach them. Note that the optimized inner radius 
of the ring-shaped beam is generally found to be of the order of the main beam size Wq (so 
that the cold parts of the mirror are effectively heated by the compensation beam). 



VIII. THE CASE OF THE ADVANCED LIGO CONFIGURATION 

It is interesting to state the potential of compensation methods with a 2nd generation 
configuration for which the design is well advanced and precise numbers are available. For 
example, the "advanced LIGO" design plans either Silica substrates with absorption about 
0.5 ppm/cm, mirror dimensions a ~ 19.4 cm and h ^ 15.4 cm, a beam size Wq ~ 5.5 cm, 1.4 
kW crossing the substrates and 530 kW stored in the cavities or Sapphire substrates with 
absorption around 40 ppm/cm, mirror dimensions a ~ 15.7 cm and /i ~ 13 cm, a beam size 

~ 6 cm, 2.1 kW crossing the substrates and 830 kW stored in the cavities. For both cases, 
the coating absorption is expected at the level of e ~ 0.5 ppm, hence a factor 2 better than 
today's best coatings. With such numbers, and without thermal compensation, the losses 
for Silica substrates are about 0.15, while they are about 1.8 x 10^ for Sapphire substrates. 
It is clear again that thermal compensation is needed for this precise configuration, whatever 
the substrate materials, unless absorptions can be drastically lowered. It is not likely to be 
much the case for Silica, but prospects for Sapphire are more optimistic, and absorption 



around 20 ppm/cm or less, as already observed in small samples |]19[, is not out of reach. 
With thermal compensation using an extra ring-shaped beam, the minimal losses achievable 
with Silica substrates are found to be about 5 x 10^^ (less than 2 x 10^^ with removal 
of the parobolic part). This result means already (at least) losses around 1% only due to 
thermal lensing in input mirrors of the kilometric cavities. This is not encouraging for using 
Silica substrates. With Sapphire substrates now, the figure 12 shows the minimal losses as a 
function of the absorption in the range 10-50 ppm/cm (since it is today like a free parameter 
that will probably evolve favourably). Already with absorption about 40 ppm/cm. Sapphire 
does better than Silica, with losses about 10^^ (2 x 10^'^ with parabolic part removal). With 
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absorption about 20 ppm/cm, the losses decrease down to about 460 ppm (some tens of ppm 
with parabolic part removal). With Sapphire substrates and thermal compensation, thermal 
lensing should not be longer a source of concern for the good working of the interferometer. 

IX. CONCLUSIONS 

In this paper, we have theoretically investigated several means to compensate for thermal 
effects in mirrors of (advanced) interferometric detectors of gravitational waves. The best 
method seems to heat the cold parts of the substrates by absorption of a ring-shaped beam, 
while the method of thermostating the circumference of the mirror is clearly not competitive. 
With the VIRGO numbers. Silica substrates and absorption coefficients e = 1 ppm and a = 1 
ppm/cm, and in case of increase of the circulating powers by a factor 10, the losses can be 
lower down to £ ~ 2.7 x lO"^, vs £ ~ 0.16 if nothing is done. 

A factor about 2 can still be gained if we can remove the parabolic part of the thermal 
lens, but this supposes that the thermal lenses are well symmetric between the two arms 
of the interferometer. For a beam size wo — 2 cm, the order of magnitude of the light 
source power needed to achieve the minimal losses value is larger than about 40 W, since an 
absorbed power of about 38 W is required. This needed power can be lowered with a larger 
beam size. For instance, only 4.5 W are needed to be absorbed ii Wq — 6 cm, that means 
a gain of almost one order of magnitude for the required power. Of course the beam size 
cannot be increased much further due to the finished size of the mirrors (typically not larger 
than 20 cm). 

Now, again with the example of the VIRGO numbers, if we want a gain of one order of 
magnitude for the shot-noise limited sensitivity, that means an increase of the circulating 
powers by a factor 100 this time. If the absorption coefficients in the coatings and in the 
Silica substrates can not be decreased by at least one order of magnitude, it is then clear 
that we will have to give up the Silica for the substrates, in favour of Sapphire for instance, 
since it will no longer be possible to compensate for thermal gradients (the compensation 
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methods can not do miracles: eg with Pcoat ~ 1 MW, Pgub ~ 100 kW, e—1 ppm and a — 1 
ppm/cm, the losses can not be lower than 97% !). This means finally that the mirror R&D 
effort needed for future upgrades has two options: either going on with Silica susbtrates and 
improve absorption coefficients at 1.064 //m (better coating technology, lesser OH radicals in 
the substrate ...) or stopping with Silica technology and developping substrates with other 
materials like Sapphire. Numerical results obtained with the advanced LIGO configuration 
are unambiguous regarding the substrate materials : Sapphire is (will) be better. 
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FIG. 1. Heating of a cylindrical mirror by a Gaussian beam. The mirror can absorb light power 
in its coating (cavity side) and in its substrate. 
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FIG. 2. Upper left: temperature field for 1 W absorbed in the coating. Upper right: temperature 
field for 1 W absorbed in the substrate. Lower left: the thermal lens profile (optical path difference) 
for 1 W absorbed in the coating. Lower right: the thermal lens profile for 1 W absorbed in the 
substrate. 
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FIG. 3. Upper left: temperature field for 15 mW absorbed in the coating and 10 mW in the 
substrate (typical VIRGO numbers). The temperature scale is in K. Upper right: temperature 
field for 150 mW absorbed in the coating and 100 mW in the substrate (VIRGO upgrade numbers). 
The temperature scale is in K. Lower left: the thermal lens profile (optical path difference) for 15 
mW absorbed in the coating and 10 mW in the substrate (typical VIRGO numbers). The OPD 
scale is in /xm. Lower right: the thermal lens profile for 150 mW absorbed in the coating and 100 
mW in the substrate (VIRGO upgrade numbers). The OPD scale is in /xm. 
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FIG. 5. Losses a as a function of the power absorbed from the extra beam. The losses are 
minimal for the quite large power Pabs — 38 W. 
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FIG. 6. Thermal lens profiles for respectively Pabs — 3.8VF (10 times too low), Pabs — 38W 
(optimal case) and Pabs — 380VF (10 times too high). In the optimal case, we note the flatness of 
the lens in the central part of the mirror. 
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FIG. 7. Losses a as a function of the power absorbed from the extra beam. The parameters 
are ai = 2 cm, 02 = 20 cm and 0.2 = 0.2/cm. The losses are minimal for the quite large power 
Pabs ^ 38 W. 
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FIG. 8. RMS of the distribution of X] 2Jo(/cm'")/CmJi(Cm) for a sampling of r in [0, a] as a 
function of iV, number of terms in the series. The RMS is minimal for N 2± 700. 
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FIG. 9. Losses a as a function of the thermostat temperature. The tower temperature is fixed 
to Te ~ 300 K. 
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FIG. 10. Top: temperature field in the mirror. Bottom: thermal lens profile. Parameters: 
Te = 300 K and Tq = 311 K (optimum case). 
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FIG. 12. Losses (in ppm) vs the Sapphire substrate absorption. For each absorption, the losses 
are optimized with respect to the inner radius of the compensation beam and to the absorbed 
power from this beam. 
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